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ABSTRACT 

We present a study of accretion in a sample of 45 young, low mass objects in 
a variety of star forming regions and young associations, about half of which 
are likely substellar. Based primarily on the presence of broad, asymmetric 
Ha emission, we have identified 13 objects (~30% of our sample) which are strong 
candidates for ongoing accretion. At least 3 of these are substellar. We do not de- 
tect significant continuum veiling in most of the accretors with late spectral types 
(M5-M7). Accretion shock models show that lack of measurable veiling allows 
us to place an upper limit to the mass accretion rates of < 10~^° Mq yr~^. Using 
magnetospheric accretion models with appropriate (sub)stellar parameters, we 
can successfully explain the accretor Ha emission line profiles, and derive quan- 
titative estimates of accretion rates in the range 10~^^ < M < 10~^ MqJt~^. 
There is a clear trend of decreasing accretion rate with stellar mass, with mean 
accretion rates declining by 3-4 orders of magnitude over ~ 1 — 0.05 Mq. 

Subject headings: accretion disks, brown dwarfs, stars: emission-line, pre-main 
sequence, circumstellar matter, 

1. Introduction 

The collapse (~ 10^ yr; Adams, Lada, & Shu, 1987) of a slowly rotating molecular cloud 
core leads to formation of a protostar and surrounding circumstellar accretion disk which can 
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last up to ~ 10^ yr in at least some cases (MuzeroUe et al., 2000a; Alencar & Batalha, 2002). 
Observable signatures of disks and disk accretion include far- and mid-infrared emission from 
dust at a range of temperatures, near-infrared emission from hot dust and gas in the inner 
disk, and a variety of emission lines and UV/optical continuum excess due to accretion of 
material from the inner disk directly onto the star (see, e.g. Mannings, Boss, & Russell 2000 
for reviews). 

Historically, accretion from the disk onto the star was thought to occur through a hot 
equatorial "boundary layer" (e.g. Lyndcn-Bell & Pringlc 1974; Bertout, Basri, & Bouvier 
1988; Basri & Bertout 1989). Over the past decade, however, a new model of magneto spheric 
accretion has been widely adopted for classical T Tauri stars (hereafter, CTTSs). This model 
was originally discussed in the context of neutron stars by Ghosh & Lamb (1979) and first 
applied to CTTSs by Uchida & Shibata (1985), Bertout et al. (1988), and Konigl (1991). 
In this scenario, the inner regions of the circumstellar disk are disrupted by the stellar 
magnetosphere, which channels viscously accreting material out of the disk plane and towards 
the star along magnetic field lines; an accretion shock is formed as the free-falling material 
impacts the stellar surface. This model can also generate magnetocentrifugal outflows, which 
are launched from a region just exterior to the point of disk truncation (Shu et al, 1994), or 
from near the poles of the star (Hirose et al. 1997). Several of the detailed observational 
characteristics of CTTSs, such as blue continuum excesses including the Paschen and Balmer 
continuum slopes (e.g. Calvet & GuUbring, 1998), relatively small near- infrared excesses 
(e.g. Meyer, Calvet & Hillenbrand, 1997), permitted emission-hne profiles (e.g. MuzeroUe, 
Calvet, & Hartmann, 2001a; hereafter, MCH), magnetic field strengths ( Johns- KruU et al. 
1999; Johns-KruU & Valenti 2000), and slow rotation periods (e.g. Edwards et al. 1993), 
are successfully explained by such a model. 

The presence of inverse P Cygni hne profiles in the spectra of many CTTSs motivated 
Calvet & Hartmann (1992) and then Hartmann, Hewett, & Calvet (1994) to first investi- 
gate magnetospheric accretion flows as the source of the line emission. Radiative transfer 
calculations of magnetospheric accretion predict a characteristic blue asymmetry - with or 
without a redshifted absorption feature - in both hydrogen (Balmer, Paschen, and Brackett 
lines) and metallic (notably NaD, 01, and the Call triplet) line proflles. These calculations 
successfully reproduce both the line fluxes and the line proflles (Hartmann, Hewett, & Cal- 
vet 1994; MuzeroUe et al. 1998ab; MCH) observed in most CTTSs in the Taurus molecular 
cloud. 

Results thus far, however, have been limited to T Tauri stars of restricted range in 
mass and age (M ~ 0.5 M©; r <l-3 Myr ). Few high-resolution observations even exist 
of stars with masses < 0.2 Mq. In order to develop a broader understanding of accretion 
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disk activity, in this contribution we extend the analysis of accretion to a mass range more 
representative of the stellar mass function, focusing on young brown dwarfs and stars in the 
mass range 0.05-0.2 Mq (hereafter referred to as very low mass, or "VLM", objects). Our 
results for a single object among this sample, V410 Anon 13, appeared in Muzerolle et al. 
(2000b), the first measurement of accretion in a candidate sub-stellar mass young T Tauri 
object. The derived M of ~ 5 x 10~^^ M0yr~^ remains the lowest measured accretion rate, 
even when the larger sample presented here is considered. We have also recently extended 
the analysis in age to study ~10 Myr old stars via an investigation of disk accretion in the 
TW Hya association (Muzerolle et al 2000a; Muzerolle et al. 2001b). 

We describe our sample in section 2 and our observations, moderate- to high-resolution 
optical spectra, in section 3. Section 4 presents results from the spectra, including Ha line 
profiles, measurements of continuum veiling, and model-derived mass accretion rates, as well 
as radial and rotation velocities for part of the sample. Finally, we discuss the implications 
of our results for accretion in young VLM objects in section 5. 



2. Sample 

The sample was selected from previously identified young low mass objects in several 
different star forming regions and young clusters, including Taurus, IC 348, a Ori, Upper 
Sco, and p Oph (see references in Table 1). Our observations were designed in early 2000 
to include every known object with spectral type M5 or later and age ~10 Myr or younger. 
Since our observations were obtained, many more such objects have been discovered in star 
forming regions and young clusters (e.g. Briceno et al. 2002) and thus our sample is not 
complete, though it is likely to be representative of these very low masses. The final observed 
sample encompasses 45 young stars and brown dwarfs with spectral types M4-M8.5. At the 
0.1-8 Myr ages of these regions, the masses inferred from pre-main sequence evolutionary 
tracks (e.g. D'Antona & Mazzitelli, 1998; Chabrier et al. 2000) span the substellar boundary, 
ranging from 0.02-0.14 M0. 

Table 1 contains basic information on our sample, including spectral types found in 
the literature from which, in combination with literature photometry, the other parameters 
(logT, logL, extinction, H-K excess) have been calculated. Masses and ages were determined 
from the D'Antona & Mazzitelli (1997, including unpublished updates provided in 1998) 
pre-main sequence evolution tracks, adopting the transformations used in Hillenbrand & 
Carpenter (2000). Our observations do not include enough dwarf and giant spectral tem- 
plates to re-examine with any great accuracy the previously published spectral types. We 
note that one particular object, LkHa 358, has types from Cohen & Kuhi (1979; M5.5) 
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and Kenyon et al. (1998; M5-M6) that are clearly inconsistent with our spectrum, which 
appears to be more like K8. This star also displays spatially extended Ha and [SII] doublet 
emission, likely originating in a jet. Our spectral type estimate of another object, Haro 6-28, 
is also discrepant with the previously reported spectral type of M5 (Cohen & Kuhi 1979); 
we estimate M2.5. Both of these objects exhibit accretion signatures, including continuum 
veiling. Figure 1 shows the HR diagram for our sample, and illustrates the estimated mass 
and age ranges spanned in our study of accretion diagnostics. 

3. Observations 

We obtained spectra at Keck II with the ESI echellettc spectrograph (Sheinis et al. 
2002) on Jan. 25-28, 2001, and at Keck I with the HIRES cchelle spectrograph (Vogt et al. 
1994) on Jan. 10-11 2000 and June 30, 2001. Extraction of one dimensional spectra from the 
two-dimensional images was performed using the MAKEE software written by T. Barlow. 

ESI provides wavelength coverage from 3900 - 11000 A, but with only moderate resolu- 
tion, R ~ 8000 (37.5 km s~^), while HIRES has a smaller wavelength range, for our setup, 
6320 - 8730 A with gaps between the orders and typical echelle resolution of i? ~ 34, 000 
(8.8 km s~^). All sets of data cover the important Ha emission line, our primary diagnostic 
of accretion. Despite lower resolution, the ESI data still resolve this line when accretion is 
present, since in that case the line widths are typically several hundred km s~^; however, 
chromospheric emission line profiles and photosphcric atomic absorption features generally 
are not resolved. The trade-off for lower resolution is the benefit of increased wavelength 
coverage into the blue which yields a larger baseline for measuring continuum veiling (al- 
though in most of our objects, the S/N is too low to make measurements below ~ 5000 
A). Line profiles are more easily resolved with HIRES, but gaps between the orders prevent 
simultaneous measurement of all emission hnes of interest for accretion studies. 



4. Results 

4.1. Radial and Rotational Velocities 

We determined v sin i values for the HIRES sample only - the resolution of ESI is too 
low for any meaningful limits given typical values for low mass young stars (e.g., Hartmann 
et al. 1986; White & Basri 2003, hereafter WB03). The rotation velocities were determined 
by cross-correlating photospheric absorption features in small portions of each spectrum (see 
Tonry & Davis 1979; Hartmann et al. 1986). We selected intervals of ~ 20 A in several 
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different orders containing relatively strong absorption features with good S/N. Values of 
vsini were derived from the widths of the cross-correlation peaks, which were calibrated 
using template spectra of Gl 447 {v sin i < 2 km s~^, unresolved at HIRES resolution; Dclfosse 
et al. 1998) convolved with a range of rotation velocities. Results for the different orders 
were then averaged together, with the standard deviation taken as the final uncertainty. Wc 
present these values in Table 2. The rotation velocities range from 8-35 km s^^, with most 
being < 20 km s~^. The two fastest rotators, FP Tau and LkCa 1, both at 34 km s~^, show 
evidence of line doubling in the Ca II triplet emission lines; thus, the absorption line widths 
(which do not show obvious signs of doubling) may be due to binarity rather than rapid 
rotation. 

We determined heliocentric radial velocities through similar cross-correlation techniques 
for both the HIRES and ESI samples (see Table 2). For the HIRES spectra, we used the 
dwarf standards Gl 402 and Gl 447 (K = —2 and -31 km s~^, respectively; Delfosse et al. 
1998) as templates. For the ESI sample, we used Gl 83.1 and Gl 406 (K- = -28.567 and 
19.482 km s~^, respectively; Nidever et al. 2002). Due to the relatively low resolution of the 
ESI data (R = 8,000), the cross-correlation power is dominated by broad molecular features 
which forced us to be selective in the wavelength regime and range employed. As for the 
vsini measurements, results for different wavelength ranges were averaged together, with 
the errors taken from the standard deviation. 



4.2. Membership Information 

Many of the stars we have studied with HIRES and ESI are candidate (as opposed 
to confirmed) members of the Taurus, IC 348, Sigma Ori, Sco-Cen, and Ophiuchus stellar 
associations. In these cases, membership is supported by consistency in a color-magnitude 
diagram with the expectations for a low-mass member at the distance and age of the cluster, 
combined with a late spectral type derived from low dispersion spectroscopy. Our higher 
dispersion data contain a variety of spectroscopic indicators that may further inform mem- 
bership status. At the 1-10 Myr ages of the associations contained in our sample, high 
lithium abundance and low surface gravity are expected (we list Li equivalent widths in 
Table 2). Further, we have checked consistency of measured radial velocities with those 
previously published for known cluster members. 

In Taurus, we find that CIDA-13 lacks substantial Li I, although its radial velocity is 
consistent with the average for Taurus members; it is likely a foreground dwarf star. The 
Sigma Ori cluster has two candidate members (SOri-45 and SOri-46) with K 1 lines that 
are more similar to those of M dwarfs than of low-gravity pre-main sequence stars, as well 
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as radial velocities which are substantially different from the other Sigma Ori stars. In the 
Upper Sco association, we find three stars (UScoCTIO-85, UScoCTIO-121, UScoCTIO-132) 
which both lack substantial Li I and have dwarf-like K I absorption. 

As suggested in the tablenotes to Table 1, we question the membership of these 6 stars 
in their respective associations. Many of these stars are among the weakest Ha emitters 
(Table 3), and the emission line characteristics in general are more consistent with active 
main sequence dwarfs rather than young TTSs/brown dwarfs. We do not include these 
objects in any of the subsequent statistics or analysis. 

4.3. Emission Lines 

Emission line measurements for our entire sample are listed in Table 3, including equiv- 
alent widths for a variety of permitted and forbidden lines, as well as the Ha line widths 
at 10% of peak intensity (a good discriminator of accretion versus chromospheric emission; 
Muzerolle et al. 2000a, WB03). Measurements for two main sequence late-M dwarfs (Gl 
406 and LHS 2351) arc also given as a comparison. The emission line spectra of the VLM 
young objects arc very similar to those of higher- mass young stars, including both CTTSs 
and "weak line" TTSs (WTTSs). Balmer and He I emission are virtually ubiquitous, and 
are typically stronger than in the main sequence dwarfs. The equivalent widths are larger 
than seen in higher-mass WTTSs, but this is probably a result of the cooler photospheres, 
rather than intrinsically stronger emission. Ca II triplet emission is detected in nearly half 
of the sample, while forbidden emission such as [01] A6300 is detected in only ~25% of the 
sample (though this smaller fraction may be due in part to the drop-off in S/N of many of 
the spectra around this wavelength). 

Most of the emission line profiles are narrow (< 200 km s~^) and symmetric, typical 
of WTTSs and dMe stars. However, 13 objects show broader profiles at Ha (see Fig. 2), 
including many with blueshifted asymmetries, and several with blueshifted absorption com- 
ponents (e.g., FP Tau, MHO-6, IC348-415). Such features are more reminiscent of CTTSs. 
The shapes of the upper Balmer lines generally mirror Ha , while somewhat narrower, and 
blueshifted absorption components decrease in strength as one moves up the series. Broad 
emission components in the Ca II and He I lines are seen only in the earliest-type object in 
the sample, LkHa 358. This star is also the only one that exhibits a high- velocity component 
at [01]. Overall, the most active of the 0.05 — 0.2 Mq VLM young objects have emission 
characteristics similar to the least active among the well-studied 0.5 Mq CTTSs. 
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4.4. Accretion Signatures 

We interpret our data in the context of accretion disks and accretion-driven outflows 
associated with these young VLM objects, following our previous work on higher mass TTS. 
Hot UV/optical continuum veiling of photospheric absorption features, as well as narrow 
emission components in Ca II and He I, provide clues to the density, temperature structure 
and filling factor of the accretion shock; these parameters may differ from those characteristic 
of the well-studied 0.5 M© stars, as M*/ varies. Observations of broad Balmer lines provide 
detailed information on the nature of accreting magnetospheres. Forbidden emission lines 
such as [01] and [SII], and blueshifted absorption components in permitted lines, are typically 
associated with accretion disk-driven outflows (e.g. Hartigan, Edwards, & Ghandour 1995; 
Calvet 1997). We use all of these features to identify objects which are likely undergoing 
active accretion, but focus primarily on veiling and Ha emission line profiles in our analysis. 

Ji..4-1- Veiling 

One familiar manifestation of accretion in T Tauri stars is the veiling of photospheric 
absorption features by excess continuum emission produced in the accretion shock. Veiling is 
typically defined as the ratio of the excess flux to the photospheric flux at a given wavelength: 
rx = Fx^ex/ Fx^cont- Consistent with the historical suggestions that accretion signatures are 
weaker in lower mass stars, we find very little veiling in our VLM sample. We have measured 
or placed upper limits on the veiling for each of the 13 likely accretors (as defined from 
emission line indicators; see below), using comparison spectra of similar spectral type (e.g., 
Hartigan et al. 1989). In order to avoid mismatches due to gravity-sensitive features (see 
WB03 for a more detailed discussion), we opted to use the non-accreting subsample, defined 
as those which clearly lack infalling gas signatures at Ha , as veiling templates. Since veiling 
at red optical wavelengths is much less sensitive a diagnostic of accretion than Ha , these 
templates should not have any intrinsic veiling of their own. We note that WB03 measured 
an upper limit r < 0.22 at 6500 A for one of our templates, V410 X-ray 3, using a combination 
of dwarf and giant standard spectra. 

Veiling measurements at several wavelength regions are given in Tables 4 and 5. Due to 
the gaps in the wavelength coverage of the HIRES spectra, we could not measure veiling at 
exactly the same wavelength regions as ESI, hence we list these separately. The wavelength 
regions chosen for the ESI sample contain broad molecular absorption features; we focused 
on these rather than narrow atomic features in order to avoid line blending problems due 
to the lower resolution of the ESI spectrograph. For the two objects earlier than ~M3, 
Haro 6-28 and LkHa 358, the derived veiling is more uncertain since the atomic features 
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dominate. These problems are much less significant for the HIRES sample; however, the grid 
of templates is less extensive, and slight mismatches due to differing rotation rates result in 
greater uncertainties. Note especially the results for the V410 Anon 13 and CIDA 14 HIRES 
spectra - the closest template in spectral type is about 0.5 subclasses earlier, and rotates 
faster, which could result in slight underestimates of the true veiling (although the results 
for V410 Anon 13 arc consistent with those of the ESI spectrum). 

Most of the objects at M4 or later exhibit no significant veiling (i.e., > 2a of the typical 
uncertainty of ±0.1) at any of the wavelengths. One low-mass object, IC348-415 does have 
significant veiling with r ~ 0.5 at 6200 A, though this value is highly uncertain due to low 
S/N; this object also has the strongest Ha emission of the entire sample, and our best profile 
model indicates the highest M of the >M4 sample. Nevertheless, we consider this detection 
to be marginal at best. Only in the 2 earhest-type objects in the sample, Haro 6-28 and 
LkHcc 358, do we measure robust levels of veiling, indicating higher levels of accretion than 
in the lowest-mass sample as a whole. The veiling levels are consistent with the general 
emission line characteristics of the stars, such as the presence or absence of broad Ca II 
triplet emission components, which were noted by Muzerolle, Hartmann, & Calvet (1998) 
to have a threshold behavior whereby the broad component does not appear when the mass 
accretion rate is < 10"® M0yr~^. 



44.2. Ha Profiles 

Continuum-normalized Ha line profiles are shown in Figure 2. Traditionally, the thresh- 
old between CTTSs and WTTSs has been set at a canonical boundary of Ha EW = 10 A, 
which divides accreting from non-accreting TTSs around the peak of the mass function in 
Taurus, ~ 0.5 Mq. However, this threshold is inappropriate for much lower mass objects, 
where the decreased optical photospheric continuum results in typically larger equivalent 
widths despite comparable emission fluxes from chromospheric activity. Recently, WB03 
compared Ha 10% line widths to equivalent widths for TTSs of a wide range of masses, and 
calibrated a spectral type-dependent equivalent width boundary between chromospheric and 
accretion emission at Ha , down to late-M objects. Figure 3 shows the distribution of Ha 10% 
line widths, as well as a comparison with equivalent widths, for our bona- fide young object 
sample. Here, we also focus on the fact that the line profile shape is a much more robust dis- 
criminator of chromospheric vs. accretion emission. The broad (>200 km s~^), asymmetric 
profiles produced in infalling accretion fiows are easily identified even in the lower-resolution 
ESI spectra. Ha emission from WTTSs, on the other hand, is similar to that seen in chro- 
mospherically active main-sequence dwarfs (e.g., Worden, Schneeberger, & Giampapa 1981), 
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with relatively small equivalent widths and narrow (<200 km s~^), symmetric line profiles. 

Our accretion models (see below) indicate that accretion profiles in VLM objects should 
be detectable down to very low accretion rates (~ 10~^^ M0yr~^), much lower than the 
sensitivity limits of any other optical diagnostics, including continuum veiling^. 

We have identified hkely accretors on the basis of line profile shape, namely whether 
the Ha lines are broad (> 200 km s~^) and asymmetric. A total of 10 objects meet both of 
these criteria, and we consider them to be strong candidates for accretion. Another 4 objects 
(V927 Tau, IC348-205, 1C348-382, and UScoCTIO-75) meet the line width criterion, but are 
relatively symmetric. Yet another 3 objects (MHO-4, MHO-5, and UScoCTIO-128) meet 
neither of the Ha accretion criteria, but show other potential signposts of accretion, such as 
forbidden line emission or Ca H triplet emission fluxes stronger than typical dwarf /WTTS 
levels. Of these 7 objects with ambiguous accretion indicators, we beheve at least 3 are 
compeUing cases for the presence of active accretion. IC348-205 and IC348-382 have consid- 
erably stronger Balmer and Ca H triplet emission than the VLM WTTSs, and, in both cases. 
Ha is quite broad, if relatively symmetric, with two of the largest EWs in the entire sam- 
ple. MHO-5 exhibits relatively narrow and symmetric Ha emission; however, the equivalent 
width is larger than typical of VLM WTTSs of similar spectral type (e.g., WB03). The pro- 
file shape of Ha in this object is probably due to a pole-on geometry for the accretion flow, 
which produces narrower, more symmetric emission profiles since most of the high-velocity 
infalling gas is moving perpendicular to the hue of sight (§5.2; MCH). MHO-5 also exhibits 
surprisingly strong [01] emission, which is never seen in higher-mass WTTSs. Briceiio et al. 
1998 detected similarly strong Ha and [01] emission in MHO-5, so we discount variability 
due to flaring. For our subsequent analysis, we consider these 3 objects, plus the original 
10 Ha -selected strong candidates, as our total "accretor" sample (see Tables 4 and 5). The 
distribution of Ha 10% line widths and equivalent widths for this subsample are delineated 
in Figure 3. The remaining 22 objects exhibit emission characteristics consistent with chro- 
mospheric activity alone - narrow, symmetric Ha , no forbidden emission, and little or no 
Ca II emission. 

Yet another ambiguous case of interest is V927 Tau, whose Ha profile includes a narrow. 



^There is a possibility that at least some of the broad Ha profiles could be due to flares. We believe 
this is unlikely for several reasons. Previously published profiles of M-dwarf flares sometimes show a very 
weak broad gaussian component {FWHM ^ 200 km s~^), often redshifted, superposed on the (stronger) 
chromospheric profile, and the equivalent widths are typically < 20A (e.g., Abdul- Aziz ct al. 1995; Abranin 
et al. 1998; Martin 1999). V927 Tau, which shows a similar profile, could be such a case. In general, 
however, these flare proflles are not consistent with our observations. Also, our repeat observations of V410 
Anon 13, separated by 9 months, show a nearly identical Ha profile, indicating a relatively stable source for 
the emission, thus favoring accretion over flare activity. 
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symmetric component with a small self-absorption reversal at line center, all hallmarks of 
chromospheric emission. However, there is also a very weak, slightly asymmetric broad 
component that could be due to accretion. The equivalent width of the broad component is 
roue; hly 2 A, more than an order of magnitude lower than is typical of the accretor sample. 
These characteristics are reminiscent of flare activity, and we cannot reliably distinguish 
between that or very low-level accretion without repeat observations, so we do not include 
this object in the accretor sample. 



5. Modeling Accretion Flows 

5.1. Continuum Excess Models 

Following the procedures of Calvet & GuUbring (1998, CG98), we calculate accretion 
shock emission and estimate the amount of veiling expected for a given mass accretion rate. 
In short, inf ailing magnetospheric material accretes onto the stellar photosphere through a 
shock just above the photosphere. The material is assumed to be in plane-parallel columns. 
In each column, soft X-ray emission from the shock heats both the pre-shock region above 
the shock and the stellar photosphere below it. The shock flux Fghock is essentially the sum 
of the spectra from the heated photosphere and the preshock region. The spectral shape of 
the shock flux depends on the energy flux carried by each accreting column, JF, while the 
level of the emission depends on the surface coverage of accreting columns, f^TiRl (CG98). 
The total luminosity of the shock is 

Lshock = f^TTRlT, (1) 

neglecting the intrinsic stellar flux, which is at least an order of magnitude lower than the 
flux carried by the column in these cool stars. At the same time, 

where — 1 — 1/Rmag, with Rmag being the radius at which magnetospheric infall begins. 

We have calculated models for two generic stars with parameters corresponding to our 
sample: an M5 star (Te// = 2900 K, = 0.12 Mg, = 1.2 Rq) and an M6 star {T^ff = 
2800 K, = 0.05 Mo, R^ = 0.5 Rq). Models were calculated for log T = 10, 11, and 
12, which correspond to the range of energy fluxes that provide better fits to the optical 
and ultraviolet excess fiuxes in higher mass T Tauri stars (CG98; Gullbring et al. 2000). 
The emission from the shock was calculated for M — 10~^^ M0yr~^,lO~^ Moyr"^ and 
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10 ^ Mq yr ^, using eqs. (1) and (2) to obtain the filling factors / for given values of M and 
We take Rmag = 3, in accord with the Ha emission line models (§6). 

To calculate the veiling, we compare the flux arising from the accretion shock with the 
intrinsic photospheric flux. Here, we calculated the photospheric flux using bolometric cor- 
rections from Luhman (1999) to convert from model luminosities to absolute J magnitudes, 
and then used standard colors to obtain corresponding V, Rc, and Ic magnitudes. Colors 
were taken from Lcggctt (1992) for types later than M4, and from Kenyon & Hartmann 
(1995) for earlier types. Fluxes at the cfTcctivc wavelength of each filter were obtained using 
the zero points of Bcsscll & Brett (1988). Fluxes were then scaled to the assumed model 
radii, and then the veiling was calculated at each filter wavelength by 

^ Fexcess f F shock 

Pphotosphere (1 f^Fphot 

In Figure 4, we show the veihng expected for the two model stars for log M — -10, 
-9, and -8, and log T — 10, 11, and 12 (note the logarithmic scale). For a given M, the 
veihng at 5500 A does not change much with though the model with lower T has a larger 
veiling at longer wavelengths. This arises from the fact that at optical wavelengths, the 
shock flux is dominated by quasi-blackbody emission from the heated photosphere, which 
has a temperature T^p ~ [2>J- / Aa)^/^ (CG98). Thus, models with lower T emit more at 
redder wavelengths. 

It is apparent from Figure 4 that the negligible veiling measurements for most of the 
accretors (Tables 4 and 5) indicate M < 10~^° Mq yr^^. In Figure 5, we compare the detected 
veiling for 4 accretors with predictions from shock models. Energy fluxes that provide the 
best fit to the observations, as well as the estimated mass accretion rate are also shown. 
The filling factors of the shock, given in the caption of Figure 5, are rather small, << 1%, 
similar to the case of the more massive CTTSs (CG98). The accretion rates are in excellent 
agreement with those determined from Ha models, as described in the next section. 



5.2. Ha Emission Models 

Previously, we have developed radiative transfer models of Balmer emission line profiles 
for ~ 0.5 Mq CTTS (MCH). These models treat line emission produced by magnetospheric 
accretion flows, where gas accreting through a circumstellar disk is channeled onto the star 
by the stellar magnetic held. The model assumes an ideal dipolar field geometry, and a 
specified gas temperature distribution and mass accretion rate; the gas density distribution 
is determined directly from the geometry and the gas velocity (a function of the mass and 
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radius of the star) . Further details can be found in MCH. Typical results show a characteristic 
profile shape with blue asymmetry, due to occultation effects and some absorption of the 
red wing where the receding part of the flow is seen in projection against the star and 
accretion shock. Frequently, outright redshifted absorption components are seen in the upper 
Balmer lines, but rarely in Ha because of thermalization and opacity broadening effects. The 
qualitative similarity of most of the asymmetric Ha profiles of the sample presented herein to 
these general model characteristics lends strong support for our application of these models 
to the lower-mass objects. 

Here, wc employ the same model, using two values of the stellar mass (0.05 and 0.15 
Mq) and radius (0.5 and 1 Rq) that span the range appropriate for our VLM sample (see 
Table 1). The mass range will change somewhat depending on which evolutionary tracks are 
used; however, this should not affect our modeling applications since the model gas velocities, 
and, hence, line widths, are only sensitive to y/M, and the derived masses are unlikely to be 
off by more than a factor of two. The radius is dominated by uncertainties in the spectral 
type - T^,ff transformation; a 200 K error in Te// translates to ~ 30% uncertainty in the 
radius. Since the model gas density does not depend very strongly on the radius (see MCH, 
and references therein) , our model constraints are insensitive to this level of uncertainty. 

We use the same temperature constraints as those found to reproduce the 0.5 Mq CTTS 
profiles, which should be a reasonable assumption since the gas heating is likely mechanical 
and independent of the stellar T^^ff. These constraints result in an inverse proportionality 
between the gas temperature and density. Since the accretion rates of all the objects wc arc 
considering here are < 10"^ MQyr"^ (as inferred from the lack of continuum veiling), the 
accretion flow gas temperatures must be > 10,000 K to obtain sufficient Ha emission, thus 
the gas is almost completely ionized. As an independent check of the adopted temperature 
constraints, we show several Ha models as a function of both temperature and M in Figure 6. 
The models at the higher value of M (10~^ Mq yr~^) exhibit broader line widths than those 
at the lower M, as opacity broadening effects kick in due to the higher density; the profile is 
still broader at a lower temperature of Tmax = 8000 K (temperatures lower than this are not 
applicable at this value of M as the line flux becomes much lower than observed). At the 
lower value of M (10~^° M0yr~^), the density is too low to produce opacity broadening in 
the wings, and the proflle shape is nearly identical for T^^x = 10, 000 and 12,000 K. Thus, 
the model temperature can be constrained enough to obtain more accurate constraints on 
the mass accretion rate. 

As for the flow geometry, we have selected the flducial magnetosphere inner and outer 
radius values from MCH (2.2-3 R*). A larger magnetosphere would produce higher fine 
fluxes and more strongly centrally-peaked proflles for a given accretion rate, while a smaller 
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magnetosphere would produce redward asymmetries in the line profiles due to stellar occul- 

tation effects. The fiducial size is also consistent with inner hole sizes inferred from observed 
{H - K) and {K - L) excesses (Table 1; Meyer et al. 1997; Liu et al. 2002). Varying the 
magnetospheric size by roughly 20% may still produce an adequate match to the observed 
profiles, and still be consistent with infrared excesses. This results in an uncertainty in our 
M constraints of about a factor of 2-3 (see MCH for a detailed comparison of magnetospheric 
size, accretion rate, and line flux). 

Finally, we have not included the effects of rotation in the models presented here. MCH 
showed that for Vrot ^ 20 km s^-*^, consistent with most of our f sini measurements (§4.4), 
rotation does not have a significant effect on the line profile. A few of the profiles we model 
below show evidence for faster rotation, but since we cannot measure v sini, we have elected 
not to include rotation effects for these. M constraints are not affected by this omission, 
since the model gas density is virtually unchanged for the plausible range of rotation rates 
in the low mass objects. 

For each accretor Ha profile, we constructed a model by adopting the above parameter 
constraints, varying only the mass accretion rate and inclination angle to find the best 
match to the observed profile shape and fiux. A range of accretion rates were explored, 
10""^^ < Af < 10~^ M ^yr^^, with the lower limit set where accretion emission is produced 
at an undetectable level, and the upper limit set by the veiling constraints. The inclination 
angle was varied from pole-on to edge-on orientations of the star/disk system. Our best 
model matches arc shown in Figure 7, with parameters listed in Table 6. We have not 
modeled the 3 accretors earlier than M5, due to the temperature parameter degeneracy 
problem mentioned above (we did derive accretion rates from the veiling for two of these, 
Haro 6-28 and LkHcc 358, with M > 10~^ Mgyr"^ in both cases; see Figure 5). We also do 
not show a model for V410 Anon 13, the subject of an earfier paper (MuzeroUe et al. 2000b). 
Note that in two objects shown in Figure 7, MHO-6 and IC348-415, the Ha profiles show 
blueshifted absorption components. Such features are produced in accretion-powered winds, 
which are not treated in the model calculations. In these cases, we have tried to match only 
the red and high-velocity blue wings. 

The overall agreement between the observed and model profiles is remarkably good, 
especially considering the idealized model geometry. In reality, the accretion fiows are prob- 
ably not axisymmetric, and in many cases may be broken up into discrete "streams" (e.g., 
Muzerolle et al. 2000b). The discrepancies near line center for some of the IC 348 objects 
(higher predicted fiux than observed) may be a result of neglecting rotation in the model 
calculations; we do not have rotation velocities for those objects, which could conceivably be 
rotating much faster than typical of the Taurus objects in the HIRES subsample. Previous 



-14- 



models we have generated for higher-mass stars with rapid rotation generally produce some- 
what broader, less-strongly peaked profiles with doubling at line center. In addition, the 
assumption of a large velocity gradient implicit in our use of the Sobolev approximation in 
the model calculations for the line source function break down near line center, which mainly 
probes the outermost region of the magnetospheric flow where velocities are relatively small 
and densities are correspondingly large. Thus, we may be overestimating the amount of flux 
that can escape from the outer regions of the magnetosphere. Nevertheless, these models 
provide quantitative estimates of mass accretion rates (which we show in Tables 4 and 5) 
that are accurate to within a factor of ~3-5, given the preceding caveats. Note that these 
Ha -derived accretion rates are consistent with the values and upper limits derived from the 
veiling, as shown in the previous section. Although our M estimates arc not as accurate as 
UV excess measurements, they are likely the best achievable constraints for VLM objects at 
present, given the generally small accretion rates, which yield only a small amount of excess 
emission in the UV. 



6. Discussion 

The clear signatures of accretion exhibited in roughly 30% of our late-type sample 
indirectly reveal the presence of circumstellar accretion disks, and indicate that such disks 
are relatively common in young VLM objects. Observations of infrared excesses yield further 
evidence of disk accretion, and several recent studies along these lines have shown that the 
fraction of VLM young stars and brown dwarfs with disks is indeed - and not surprisingly - 
similar to that of the higher mass TTSs (e.g., Muench et al. 2001; Liu, Najita & Tokunaga 
2002). Examination of Table 1 shows that a fairly large fraction of our sample exhibits near- 
infrared excess, indicative of circumstellar disks. The level of excess at these wavelengths 
suggests circumstellar disks with fairly small inner holes, relatively consistent with the inner 
magnetospheric radius of our Ha models. However, most of the excesses are fairly small 
(A(iJ — K) < 0.2), and we cannot deflnitively rule out spots, spectral type mismatches, or 
cooler, unresolved companions as possible sources of apparent excess. 

In any case, the general agreement between the observed and model Ha proflles we 
present here shows that magnetically mediated accretion from a disk is applicable to VLM 
young objects. Thus, the disk accretion paradigm, which is a common (if not universal) 
aspect of early stellar evolution, appears to work over a very large range of masses, from 
substellar objects with M < O.OSM© to at least the highest-mass TTS at M ~ 2 M0 (and 
possibly to even higher masses, including the Herbig Ae/Be stars, some of which exhibit 
magnetospheric infall signatures; MuzeroUe et al. in preparation). A key parameter in the 
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physical description of accretion is the mass accretion rate, which guides both the growth 
of a young star/brown dwarf and the evolution of its disk. Mass accretion rates have been 
measured extensively at typical TTS masses of 0.5 — 1.5 Mq (e.g. Valenti et al. 1993; 
Hartigan, Edwards, & Ghandour 1995; Gullbring et al. 1998; Hartmann ct al. 1998; White 
& Ghez 2001); however, in most cases the mass ranges probed have not been sufficiently 
broad or complete to effectively study any mass dependences. RebuU et al. (2000; 2002) did 
estimate M for M ~ 0.25 — 1 Mq stars in Orion and NGC 2264 using [/-band excesses, and 
found a loose correlation with M ex M. However, their photometrically-selected samples are 
incomplete at average or lower values of M, and rather large extinction uncertainties result 
in a more uncertain [/-excess calibration. 

More recently, WB03 attempted to measure accretion rates from continuum veiling in a 
sample of VLM objects in Taurus (including four objects in common with our sample). Most 
of their sample lacked measurable veiling (as does ours), yielding upper limits on accretion 
rates of the order 10~^° M0yr~^. They did find two objects ^ with M < 0.2 M0, including 
one probable brown dwarf, that exhibited considerable veiling, yielding accretion rates of 
approximately 2 — 3 x 10~^ M0yr~^. Together with our previously reported accretion rate 
for VLM object V410 Anonymous 13 (Muzerolle et al. 2000b), and values for higher-mass 
Taurus TTS, WB03 also find an approximate correlation of M oc M. However, the number 
of M values at the lowest masses that are not upper limits is small. 

With Ha modeling, we can probe to accretion rates several orders of magnitude lower 
than possible from vciling/UV excesses. Figure 8 shows our derived M values as a function 
of mass, along with the WB03 results, and values for higher-mass TTS from Gullbring et 
al. (1998) and White & Ghez (2001). Note that WB03 and White & Ghez (2001) used 
the Baraffe et al. (1998) pre-main sequence evolution tracks to determine masses; we have 
scaled those masses, as well as their values of M (since they used M — LaccR*/GM^), to 
the D'Antona & Mazzitelli (1998) evolution tracks in order to be consistent with our results. 
There is a clear correlation between mass accretion rate and stellar /substellar mass, with 
an approximate proportionality of M oc M^. Such a dependence is steeper than in previous 
investigations. 

There are several caveats that must be emphasized when interpreting the above result. 
One is that the derived (sub)stellar masses are dependent on which evolutionary tracks are 
used; for example, the Baraffe et al. (1998) tracks yield systematically larger masses, by a 
factor of ~ 1.5-2. This also affects accretion rates derived from veiling or the UV excess, since 



''Both of these, GM Tau and CIDA-1, were previously identified from low-resolution spectra as continuum 
objects, with indeterminate spectral types. 
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those methods require the stellar mass and radius in order to derive M from the accretion 
luminosity. However, we checked the mass-M comparison using Baraffe et al. (1998) masses, 
and find the same proportionality. 

Another caveat is that objects accreting at levels low enough that Ha emission from the 
accretion flow (or UV emission from the accretion shock) is undetectable, would be missed 

in all surveys conducted to date. Such an effect is probably more relevant to the VLM 
objects since the minimum M at 0.1 Mq needed to produce observable Ha emission, as 
determined from our models, is within about a factor of 2 of the lowest measured value in 
our VLM sample. At 0.5 M0, the models indicate a higher minimum M (since the higher 
mass of the central object results in higher infall velocities, which in turn results in smaller 
magnetospheric gas densities for the same accretion rate), but one which is over one order 
of magnitude lower than the lowest measured value at that mass. Thus, at least at the 
higher- mass TTSs, the observed lower limit to M is probably real^. 

A flnal caveat is that VLM objects accreting at rates higher than those measured thus 
far could be missed due to either veiling completely flUing in photospheric features, or very 
high extinctions. However, the number of so-called "continuum stars" is fairly small (6 in 
Taurus, as hsted in Kenyon & Hartmann 1995), and not necessarily all are VLM objects. 

One of three Taurus continuum stars observed by WB03, GN Tau, turned out to be M2.5; the 
other three in Hartmann & Kenyon (1995) have bolomctric himinositics higher than likely 
for accreting VLM objects. Also, there is no evidence for any correlation between extinction 
and accretion activity in TTSs (not including embedded "Class 1" objects). Furthermore, 
in the extinction- limited sample of Briceiio et al. (2002), no correlation between extinction 
and (sub)stellar mass up to Ay ~ 4 was seen, thus there is no evidence that young brown 
dwarfs are selectively more extincted than higher-mass TTSs. 

The physical basis for a mass vs. M relation is not entirely clear. In the standard 
model of viscous disk accretion (e.g., Lynden-Bell & Pringle 1974), the disk mass accretion 
rate is proportional to the disk surface density times the viscosity. The currently favored 
viscosity mechanism, the Balbus-Hawley instability (e.g., Balbus & Hawley 1998), requires 
sufficient ionization of the disk gas to operate efficiently. A likely source of ionization of 
T Tauri disks is X-ray radiation produced by magnetic activity (Glassgold, Najita, & Igea 
1997; Igea & Glassgold 1999). Interestingly, recent studies have revealed a fairly strong 
correlation between X-ray luminosity and mass in young stellar objects (Feigelson et al. 



^This is also somewhat dependent on age, since older objects at the same mass have smaller radii, which 
results in somewhat higher densities at the same M. Nevertheless, all of the objects plotted in Figure 8 are 
at roughly the same age, 1-3 Myr. See Muzerolle et al. (2000a) for results at 10 Myr. 
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2003; Flaccomio et al. 2003; Mokler & Stelzer 2002; Preibisch & Zinnecker 2002), where 
Lx on average decreases by a factor of ~ 50 — 100 from 1 to 0.1 Mq. It is possible that 
the lower X-ray flux produced by less massive young stars could result in smaller accretion 
rates by way of decreased disk ionization fractions, driving the correlation we find between 
M and mass. However, it is not clear whether the difference in Lx over the mass range 
represented in Figure 8 can create enough of a change in disk ionization levels to produce 
a relation as steep as M oc M^. On the other hand, the correlation we have determined 
may be artificially steep, due to the possible selection effects described above. In any case, 
our results should provide interesting new motivation for continued theoretical work on disk 
accretion and evolution in young stars and VLM objects. 

It should be noted that the mass-M relation we find is for objects with a relatively 
narrow range of ages, ~ 1 — 3 Myr. Thus, we cannot make any conclusions about whether 
such a relation holds at younger ages, i.e., at the embedded proto(sub)stellar phase. Most 
of the mass is probably accreted during the first few x 10^ years, either in the initial cloud 
collapse, or via FU Orionis- type outbursts, which may be different from typical viscous 
accretion and hence may be unrelated to the scenario we suggest in the previous paragraph. 
The accretion rates measured here are indicative of a more quiescent, steady accretion phase, 
and do not have much impact on the subsequent mass evolution of these objects. 

We can comment on the fraction of accretors as a function of (sub) stellar age and mass. 
We have found 13 low-mass accretors with M > lO^^"^ MQyr^-*^, all of which reside in the 
young (1-3 Myr) IC 348 or Taurus associations. Conversely, we have found no low mass 
accretors in the somewhat older (3-8 Myr) Sigma Ori or Sco Cen associations. Consistent 
with previous results on the evolution of disk accretion in higher mass, ~ 0.5 Mq, stars 
(e.g. MuzeroUe et al 2000a; MuzeroUe et al. 2001b), we find here an accretion timescale 
on the order of several Myr for objects near the substellar boundary. At mass ~ 0.5 Mq 
and age ~ 2 Myr, the fraction of stars identified as accretors is about 50% (e.g., Kenyon 
& Hartmann 1995 for Taurus). Here, although we have found only about 30% of our total 
< 0.2 Mq sample to be accreting, ~50% of the Taurus sample is accreting, albeit at lower 
rates. 

We thank Russcl White. Michael Meyer, and Kevin Luhman for helpful discussions and 
advice, and Suzan Edwards for a critical reading of the manuscript. This work was supported 
in part by NASA Origins of Solar Systems grant NAG5-9670 to NC and LH. CB received 
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(FONACYT) of Venezuela. 
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Table 1. (Sub) Stellar Properties of the Observed Sample 
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log Li,„i (L -. ) 


/ 


.7 - // 


// - K 
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Av 


A(n- K) 


ago (yr) 


mass (J\. 


CIDA 13 


M3.5 


-1.47 


13.95 


0.61 


0.20 


11.85 


0.35 


-0.08 


6.94 


0.19 


CIDA 14 


M5 


-0.64 


12.24 


0.67 


0.33 


9.41 


0.34 


-0.01 


5.10 


0.12 


FN Tau 


M5 


-0.30 




0.95 


0.50 


8.25 


1.35 


0.10 


<5.00 


0.11 


FP Tau 


M2.5 


-0.49 


11.33 


0.69 


0.25 


8.97 


0.24 


0.00 


5.88 


0.22 


Haro 6-28 


M2.5 


-0.82 




1.07 


0.74 


9.27 


1.77 


0.39 


6.23 


0.23 


LkCa 1 


M4 


-0.44 


11.07 


0.77 


0.24 


8.69 


0.00 


0.01 


5.79 


0.21 


LkHa 358 


K7-M0 


0.88 


16.05 


2.09 


1.21 


9.69 


13.6 


0.19 


<5.00 


0.23 


MHO-4 


M7 


-1.14 


14.32 








0.97 




5.06 


0.06 


MHO-5 


M6 


-1.18 


13.72 


0.63 


0.48 


10.05 


0.11 


0.11 


5.63 


0.09 


MHO-6 


M4.75 


-1.16 


13.80 


0.71 


0.42 


10.63 


0.86 


0.06 


6.31 


0.13 


MHO-7 


M5.25 


-0.97 


13.18 


0.65 


0.29 


10.15 


0.40 


-0.06 


5.62 


0.12 


MHO-8 


M6 


-1.02 


13.63 


0.74 


0.38 


9.74 


0.62 


-0.02 


5.42 


0.09 


MHO-9 


M5 


-0.47 


12.95 








2.22 




<5.00 


0.11 


V410Anon 13 


M5.75 


-1.45 


16.59 


1.24 


0.71 


10.95 


3.83 


0.12 


6.51 


0.08 


V410 Xray-3 


M6 


-1.20 


14.18 


0.71 


0.39 


10.39 


0.80 


-0.02 


5.66 


0.09 


V410 Xray-5a 


M5.5 


-1.29 


15.37 


1.24 


0.61 


10.14 


2.57 


0.11 


6.33 


0.10 


V927 Tau 


M3.5 


-0.46 


11.46 


0.86 


0.31 


8.68 


0.38 


0.02 


5.17 


0.15 


IC348-165 


M5.25 


-0.98 


16.07 


0.81 


0.50 


11.77 


2.41 


0.02 


5.64 


0.12 


IC348-173 


M5.75 


-0.95 




0.72 


0.53 


11.85 


1.42 


0.09 


5.40 


0.10 


IC348-205 


M6 


-1.11 




0.69 


0.65 


12.15 


1.21 


0.21 


5.52 


0.09 


IC348-336 


M5.5 


-1.40 


17.61 


0.91 


0.63 


13.27 


3.12 


0.10 


6.42 


0.10 


IC348-363 


M8 


-2.12 


17.95 


0.64 


0.69 


13.50 


0.00 


0.24 


5.37 


0.02 


IC348-382 


M6.5 


-1.91 


18.92 


0.87 


0.76 


13.72 


2.77 


0.20 


5.69 


0.04 


IC348-407 


M7 


-2.01 


19.71 


0.88 


1.41 


13.97 


3.55 


0.78 


5.37 


0.03 


IC348-415 


M6.25 


-2.00 


18.23 


0.67 


0.12 


14.03 


1.35 


-0.34 


6.07 


0.04 



Table 1 — Continued 



object Spec. Type log L^oi (L©) I J - H H - K K Ay A{H - K) log age (yr) mass (M©) 



IC348-454 


M5.75 


_2 


19 


17 


82 





88 


0.31 


14.31 





11 


-0.04 


6.84 





05 


IC348-478 


M6.25 


-1 


92 


18 


58 


1 


13 


0.40 


14.64 


2 


27 


-0.11 


5.98 





04 


SOri-12 


M6 


-1 


51 


16 


47 





56 


0.36 


13.28 





00 


0.00 


6.31 





07 


SOri-17 


M6 


-1 


70 


16 


95 





58 


0.40 


13.79 





00 


0.04 


6.46 





06 


SOri-25 


M6.5 


-1 


75 


17 


16 





51 


0.40 


13.76 





00 


0.02 


5.62 





04 


SOri-29 


M6 


-1 


81 


17 


23 





53 


0.34 


13.96 





00 


-0.02 


6.54 





05 


SOri-40 


M7 


-2 


09 


18 


09 





54 


0.36 


14.59 





00 


-0.05 


5.39 





02 


SOri-45 


M8.5 


-2 


62 


19 


59 





65 


0.37 


15.66 





00 


-0.09 


5.98 





02 


SOri-46 


MS. 5 


-2 


71 


19 


82 





65 


0.37 


15.66 





00 


-0.09 


6.09 





02 


UScoCTIO-66 


M6 


-1 


55 


14 


85 





60 


0.38 


11.92 





00 


0.02 


6.34 





07 


UScoCTIO-75 


M6 


-1 


64 


15 


08 













00 




6.41 





06 


UScoCTIO-85 


M6 


-1 


70 


15 


23 





58 


0.39 


11.95 





00 


0.03 


6.46 





06 


UScoCTIO-100 


M7 


-1 


78 


15 


62 





62 


0.40 


11.83 





00 


-0.01 


5.31 





03 


UScoCTIO-109 


M6 


-2 


03 


16 


06 













00 




6.71 





04 


UScoCTIO-121 


M6 


-2 


19 


16 


46 





61 









00 




6.79 





04 


UScoCTIO-128 


M7 


-2 


37 


17 


09 





63 









00 




5.47 





02 


UScoCTIO-132 


M7 


-2 


59 


17 


63 





77 









00 




6.07 





02 


GY 5 


M7 


-0 


96 






1 


13 


0.66 


10.91 


5 


00 


-0.06 


<5.00 





07 


GY 37 


M6 


-1 


94 






1 


31 


0.95 


11.99 


3 


30 


0.38 


6.64 





05 


GY 141 


M8.5 


-2 


59 









76 


0.50 


13.87 





00 


0.04 


5.93 





02 



Note. — The photometry and spectroscopy used to calculate extinction values, infrared excesses, effective temperatures and 
luminosities comes from the following sources: Taurus - 2MASS, Briceno et al. (1999), Briceno et al. (1998), Luhman &: Rieke 
(1998), Kcnyon & Hartmann (1995); IC 348 - Luhman (1999), Luhman et al. (1998), Herbig (1998), Najita et al. (2000); Sigma Ori 
- 2MASS, Bejar et al. (1999) and Bejar et al. (2001); Upper Sco - 2MASS, Ardila, Martin, h Basri (2000) Ophiuchus - Wilking, 
Greene, & Meyer (1999), Luhman h Rieke (1999), Barsony et al. (1997). 
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Table 2. Velocities & Lithium 



object vsmi EW(Li A6708) 



CIDA 13 


15.9 ± 0.5 




<0.1 


/-~ i T T ~\ A -1/1 

CIDA 14 


14.0 ± 0.7 


11.8 ± 1.8 


0.51 


FN Tau 


'X A r\ 1 r\ A 

14.9 ± 0.4 


8.5 ± 0.8 


0.52 


FP Tau 


16.8 =b 2.0 


34.2" ± 0.8 


0.58 


Haro 6-28 


16.8 ± 0.9 




0.6 


LkCa 1 


9.0 ± 1.8 


ct A An 1 r\ /~> 

34.4" ± 0.6 


0.55 


LkHa 358 (ESI) 


8.9 ± 4.6 




0.4 


LkHa 358 (HIRES) 


17.9 ± 3.1 


20.4 ± 1.1 


0.58 


MHO-4 


19.1 ± 0.8 




0.5 


MHO-5 


12.3 ± 1.2 


> > > 


0.5 


MHO-6 


13.6 ± 0.7 




0.5 


MHO-7 


16.9 ± 0.8 




0.5 


MHO-8 (ESI) 


15.9 ± 0.9 




0.5 


MHO-8 (HIRES) 


15.3 ± 1.5 


16.7 ± 2.4 


0.57 


MHO-9 


13.8 ± 0.5 




0.6 


V410 Anon 13 (ESI) 


20.4 ± 0.9 




0.5 


V410 Anon 13 (HIRES) 


15.1 ± 1.3 


9.8 ± 3.0 


0.65: 


V410 Xray 3 


14.6 ± 0.9 




0.5 


V410 Xray 5 


14.7 ± 0.4 




0.6 


V927 Tau 


16.5 ± 0.6 


13.3 ± 1.3 


0.33 


IC348-165 


8.8 ± 0.6 




0.5 


IC348-173 


13.9 ± 2.4 




0.5 


IC348-205 


15.7 ± 3.4 




0.8 


IC348-336 


13.5 ± 3.0 






IC348-363 


9.5 ± 1.7 






IC348-382 


9.1 ± 0.6 






IC348-407 


7.9 ± 2.6 






IC348-415 


14.7 ± 2.5 




1.1: 


IC348-454 


10.5 ± 1.8 




0.6 


IC348-478 


10.8 ± 1.4 






SOri-12 


29.8 ± 0.7 




0.6 
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Tablc 2 — Contimicd 



object 


\r 


V sin i 


liW [Li AOiUo ) 


SOri-17 


19.66 ± 1.7 




0.8: 


SOri-25 (ESI) 


30.06 ± 2.5 




0.6 


SOri-25 (HIRES) 


29.6 ± 2.3 


9.4 ± 1.0 




SOri-29 


27.1 ± 1.6 




0.6 


SOri-40 


32.5 ± 3.3 




0.5 


SOri-45 


22.4 ± 5.3 






SOri-46 


55.8 ± 2.9 






UScoCTIO 66 


-4.4 ± 0.6 




0.6 


UScoCTIO 75 


-5.6 ± 1.1 




0.6 


UScoCTIO 85 


-24.6 ± 0.7 




<0.1 


UScoCTIO 100 


-8.9 ± 0.6 




0.6 


UScoCTIO 109 


-3.8 ± 0.7 




0.6 


UScoCTIO 121 


-38.9 ± 1.0 




<0.3 


UScoCTIO 128 


-3.0 ± 1.6 




0.5: 


UScoCTIO 132 


-8.2 ± 1.1 




<0.4: 


GY 5 


-6.3 ± 1.9 


16.8 ± 2.7 


0.5 


GY 37 








GY 141 









Note. — Velocities in km s~^, EW in A. Null values indicate 
that measurements could not be made due to low continuum 
S/N, or, in the case of vsini, could not be made due to the 
lower resolution of the ESI spectrograph. 
"Possible spectroscopic binary. 



Table 3. Emission Lines 



object Ha 

10% width 


Ha 


up 


H7 


B.S 


He I 
A5876 


[01] 
A6300 


[01] 
A6363 


Ca n 
A8498° 


Ca H 
A8542° 


Ca H 
A8662° 


accretor? 


instrument 


CIDA 13^= 


77 


-1.6 


-2.0 


-2.0 


-2.0 


-0.4 


<-0.1 


<-0.1 


<-0.1 


<-0.6 


<-0.6 


n 


ESI 


CIDA 14 


289 


-34 












<-0.1 


<-0.2 




<-0.1 


y 


HIRES 


FN Tau 


195 


-22 












<-0.2 


-0.9 




-0.7 


n 


HIRES 


FP Tau 


418 


-32 












<-0.1 


-0.1 




-0.1 


y 


HIRES 


Haro 6-28 


347 


-48 


-30 


-25 


-25 


-2.2 


-1.2 


-0.3 


-0.2 


-0.3 


-0.1 


y 


ESI 


LkCa 1 


178 


-3.9 












<-0.1 


-0.1 




<-0.1 


n 


HIRES 


LkHa 358 


502 


-63 


-22'' 






-1.9 


-15 


-4.1 


-6.3 


-6.6 


-5.1 


y 


ESI 


LkHa 358 


477 


-62 












-2.2 


-9.1 




-7.8 


y 


HIRES 


MHO-4 


116 


-43 


-56 


-35 


-25 


-4.3 


1.2: 


<-0.1 


-0.4 


-0.5 


-0.1 


n 


ESI 


MHO-5 


154 


-60 


-44 


-29 


-19 


-3.3 


-4.0 


-0.9 


-0.4 


-0.5 


-0.1 


y 


ESI ^ 


MHO-6 


309 


-25 


-13 


-11 


-8.6 


-1.7 


-1.4 


-0.3 


<-0.1 


<-0.1 


<-0.1 


y 


ESI 


MHO-7 


116 


-9.0 


-7.5 


-3.6 


-3.0 


-0.5 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


n 


ESI 


MHO-8 


115 


-18 


-17 


-11 


-5.5 


-1.4 


<-0.1 


<-0.1 


<-0.2 


<-0.1 


<-0.1 


n 


ESI 


MHO-8 


128 


-14 












<-0.1 


<-0.2 




<-0.2 


n 


HIRES 


MHO-9 


116 


-3.5 


-4.3 


-2.8 


-2.4 


-0.5 


<-0.2 


<-0.1 


<-0.1 


<-0.3 


<-0.2 


n 


ESI 


V410 Anon 13 


270 


-29 


-20'' 


-9" 




-12.8'' 


-0.8 


<-0.5 


<-0.1 


<-0.1 


<-0.1 


y 


ESI 


V410 Anon 13 


248 


-27 












<-0.4 


<-0.3 




<-0.1 


y 


HIRES 


V410 Xray-3 


116 


-20 


-23 


-18 


-10 


-2.4 


<-0.2 


<-0.1 


<-0.1 


<-0.2 


<-0.1 


n 


ESI 


V410 Xray-5a 


154 


-19 


-22 


-15** 




-2.1 


<-0.5 


<-0.2 


-0.2 


-0.3 


<-0.1 


n 


ESI 


V927 Tau 


290 


-7.0 












<-0.1 


<-0.1 




<-0.1 


n 


HIRES 


IC348-165 


347 


-54 


-28 


-24'' 




-2.4 


-0.3: 


-0.2 


-0.3: 


-0.4 


-0.4: 


y 


ESI 


IC348-173 


347 


-86 


-50 


-10'' 


-35'' 


-7.5 


-2.2: 


-0.4 


-0.8 


-1.6 


-0.7 


y 


ESI 


IC348-205 


270 


-105 


_74b 






<-1.5'' 


<-0.5 


<-0.7 


-1.1 


-1.8 


-0.9 


y 


ESI 


IC348-336 


309 


-121 


-42'' 


-20'' 




-6.4'' 


<-2.8: 


<-l.l: 


-0.2 


-0.4 


<-0.1 


y 


ESI 


IC348-363 


117 


-13 










<-21: 


<-1.0 


<-0.1 


<-0.1 


<-0.1 


n 


ESI 


IC348-382 


232 


-70 














-0.6 


-0.9 


-0.3: 


y 


ESI 


IC348-407 


155 


-24 














<-0.7 


-0.7: 


<-l.l 


n 


ESI 


IC348-415 


272 


-152 


-40'' 








<-4.0:'' 


<-1.3 


-2.8 


-4.7 


-2.6 


y 


ESI 



Table 3 — Continued 



object 


Ha 


Ha 


H/3 


H7 




He I 


[01] 


[01] 


Ca H 


Ca H 


Ca n 


accretor? 


instrument 




10% width 










A5876 


A6300 


A6363 


A8498'' 


A8542'^ 


A8662« 






IC348-454 


193 


-23 


-32'' 


-18'' 






<-3.5: 


<-3.9 


<-0.1 


<-0.1 


<-0.2 


n 


ESI 


IC348-478 


117 


-22 


-40** 








<-9.6: 


<-1.5 


<-0.1 


<-0.2 


<-0.1 


n 


ESI 


SOri-12 


116 


-9 


-9.5 


-5.5'' 






<-l.l: 


<-0.3 


<-0.1 


<-0.1 


<-0.1 


n 


ESI 


SOri-17 


126 


-4.8 


-5.0'' 








<-3.5: 


<-2.2: 


<-0.1 


<-0.1 


<-0.1 


n 


ESI 


SOri-25 


156 


-44 


-50'' 


-36'' 




-6.0 


<-1.0 


<-0.5 


-0.2: 


-0.2 


<-0.2 


n 


ESI 


SOri-25 


94: 


-36'' 














<-0.5 




<-0.25 


n 


HIRES 


SOri-29 


154 


-15 


-21 


-16'' 




<-1.5 


<-0.2 


<-0.1 


<-0.1 


<-0.1 


<-0.2 


n 


ESI 


SOri-40 


194 


-16 


-12 








<-1.0: 


<-2.3: 


<-0.1 


<-0.3 


<-0.2 


n 


ESI 


SOri-45'= 


78: 


-21'' 


-16'' 












<-0.2 


<-0.3 


<-1.4 


n 


ESI 


SOri-46'= 


194 


-14 


-9" 












<-0.4 


<-0.4 


<-0.6 


n 


ESI 


UScoCTIO-66 


116 


-6.0 


-8.8 


-8.0 


-5.5 


-1.0 


<-0.3 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


n 


1 

ESI ^ 


UScoCTIO-75 


231 


-15 


-19 


-16 


-10 


-1.3 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


n 


ESI 


UScoCTIO-85= 


154 


-7.0 


-14 


-18'' 


-5.0'' 


-0.9 


<-0.2 


<-0.2 


<-0.1 


<-0.2 


<-0.2 


n 


ESI 


UScoCTIO-100 


154 


-12 


-23 


-20 


-12 


-2.0 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


<-0.1 


n 


ESI 


UScoCTIO-109 


116 


-10 


-22 


-12 


-11" 


-2.5 


<-0.3 


<-0.1 


<-0.2 


<-0.1 


<-0.1 


n 


ESI 


UScoCTI0-121= 


116 


-7.3 


-25'' 








<-0.4: 


<-1.2 


<-0.1 


<-0.1 


<-0.2 


n 


ESI 


UScoCTIO-128 


193 


-60 


-50 


-50'' 


-30'' 


-10'' 


<-0.3 


<-0.6 


-0.3 


-0.6 


<-0.2 


n 


ESI 


UScoCTIO-132'^ 


116 


-4 










<-3.0: 


<-0.6 


<-0.2 


<-0.3 


<-0.4 


n 


ESI 


GY 5 


177 


-17 












-0.4: 


<-0.2 




<-0.1 


n 


HIRES 


GY 37 


197: 


-20'' 














<-0.4 




<-0.4 


n 


HIRES 


GY 141 


88: 


-32'' 














<-3.0 




<-17 


n 


HIRES 


Gl 406 


113 


-7.1 


-13 


-15 


-16 


-0.9^* 


<-0.1 


<-0.1 


<-0.1 


-0.1 


<-0.1 


n 


ESI 


LHS 2:ill 




-:].!) 


-11 


-9.2 


-9.1'' 


-().(!'' 


<-0.1 


<-().! 


<-0.1 


-0.2: 


<-0.2 


11 


ESI 



Note. — Line measurements are equivalent widths in A, unless otherwise noted. Ha 10% widths are in km s~^. Colons indicate uncertain 
measurements due to poor S/N, or imperfect sky subtraction at the [01] lines. Accretors have been selected based on Ha characteristics, as described 
in section 4. G1406 and LHS 2351, main sequence dwarfs with spectral types M5.5 and M7Ve, respectively, are included for comparison. 
" Ca II EW of emission component; measured against absorption pseudo-continuum when photospheric absorption component is not completely filled 
in by the emission. 

'' Continuum not detected at line position. 

Most likely not a member of the association (see text). 

Measured from a pseudo-continuum within absorption wing of Na I line. 
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Tablc 4. ESI Accrctor Sample 



object 


veiling template 


r5500 


^6200 


ttioo 


rsQOO 


logM (Moyr-i) 


Haro 6-28 


TWA 15A, CIDA 13" 


0.3 


0.2 


0.1 


0.0 


-8.7*^ 


LkHa 358 


LkCa 7 




1.0: 


1.0: 


1.0: 


-8.5'' 


MHO-5 


V410 Xray-3 


0.0 


0.1 


0.0 


0.0 


-10.8 


MHO-6 


MHO-7 


0.0 


0.1 


0.1 


0.1 


-10.3 


V410 Anonl3 


V410 Xray-3 


0.0 


0.1 


0.0 


0.1 


-11.3 


IC348-165 


MHO-8 


0.0 


0.0 


0.0 


0.0 


-10 


IC348-173 


MHO-7 


0.0 


0.1 


0.0 


0.0 


-10 


IC348-205 


MHO-8 


0.0: 


0.1 


0.0 


0.1 


-10 


IC348-336 


MHO-7 




0.2: 


0.1 


0.1 


-10 


IC348-382 


V410 Xray-3 






0.2: 


0.2: 


-10.8 


IC318-115 


VllO Xrav-3 




0.5: 


0.1 


0.1 


-9.3 



Note. — Typical veiling errors are ~ ±0.1. Colons indicate larger uncertainties due to 
poor S/N. Null values arc cases where the S/N was so low that no measurements could 
be made. Accretion rates arc in MqJT~^, all determined from Ha modeling except for '', 
which were determined from the veiling. 

" Veiling template was the average of these two stars, whose spectral types are Ml. 5 and 
M3.5, respectively. 
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Table 5. Veiling - HIRES Accretor Sample 



object 


template 


r6450 


ryioo 


r8700 


logM(Moyr-i) 


CIDA 14 


MHO-8 


-0.1 


0.1 


-0.1 


-10.3 


FP Tau 


LkCa 1 


0.0 


0.2 


-0.1 


<-9" 


LkHa 358 


LkCa 7 


0.4 


0.4 


0.4 


-8.5" 


V410 Anon 13 


MHO-8 




0.1 


0.0 


-11.3 



Note. — Typical veiling errors are ~ ±0.1. Null value indicates no 
S/N in the object spectrum around that wavelength. Accretion rates 
are in Mqjy^'^, all determined from Ha modeling except for which 
were determined from the veiling. 



Table 6. Parameters for Ha Model Comparisons 



object 


(Mq) 


{Rq) 








logM (Moyr-i) 


MHO-5 


0.05 


0.5 


2.2-3 


12,000 


45 


-10.8 


MHO-6 


0.15 


1.0 


2.2-3 


12,000 


60 


-10.3 


CIDA 14 


0.15 


1.0 


2.2-3 


12,000 


50 


-10.3 


IC348-165 


0.15 


1.0 


2.2-3 


12,000 


89 


-10.0 


IC348-173 


0.15 


1.0 


2.2-3 


12,000 


80 


-10.0 


IC348-205 


0.15 


1.0 


2.2-3 


12,000 


80 


-10.0 


IC348-336 


0.15 


1.0 


2.2-3 


12,000 


75 


-10.0 


IC348-415 


0.05 


0.5 


2.2-3 


12,000 


45 


-9.3 


IC348-382 


0.05 


0.5 


2.2-3 


12,000 


55 


-10.8 



Note. — Rmag IS thc inner and outer radius of the magnctosphcric accretion flow at 
the disk surface. T^ax is the maximum value of the temperature distribution of the flow. 
The inclination angle is measured from the axis of rotation of the system. 



-29- 




Fig. 1. — HR diagram of entire sample, using pre-main sequence tracks of D'Antona & 
Mazzitelli (1998). Solid triangles represent accretors, open triangles represent non-accretors 
(see text). 
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— Continuum-normalized Ha profiles; see Table 3 for further information. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 3. — Top panel: distribution of Ha 10% line widths for the entire sample (open his- 
togram), and for the accretor subsample (solid histogram). Bottom panel: comparison of 
Ha 10% width versus equivalent width, for accretors (solid circles) and non-accretors (open 
circles). In both plots, objects where the continuum at Ha was not detected were not 
included. 
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Fig. 4. — Logarithm of the veihng as a function of wavelength for an M5 star (upper panel) 
and an M6 star (lower panel), representative of our sample. Veiling values are shown for 
three values of the mass accretion rate, log M = -10, -9, and -8, labeled accordingly. Shock 
models are calculated for an energy flux log JF = 10 (dotted lines), 11 (sohd lines), and 12 
(dashed lines). Photospheric colors are taken from Leggett (1992) and Luhman (1999). The 
horizontal dashed lines represent the minimum measurable veiling level. 
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Fig. 5. — Veiling measurements from accretors in Tables 4 and 5 (open circles with error 
bars), compared to shock models, cf. Figure 4 (solid squares, connected). The value of 
T indicated on each panel matches the spectral shape of the veiling, while / is adjusted 
to match its level, resulting in the indicated values of M in each case (with corresponding 
errors in log M of roughly ± 0.2-0.3). The estimated filling factors / are 0.06 for LkHa 358, 
0.002 for Haro 6-28, 0.01 for IC348-336, and 0.0004 for IC348-415. 
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Fig. 6. — Model Ha profiles, as a function of magnctospheric density and temperature. All 
models calculated with = 0.05 M0, — 0.5 Rq, i — 60°, and Rmag — 2.2 — 3 R^. 
Solid line: M = lO-^M0yr-\ T^^x = 8000 K. Dashed line: M = lO-^M0yr-\ T^^x = 
10,000 K. Dot-dashed line: M = lO-^°M0yr-^ T^ax = 10,000 K. Dotted line: M = 
lO-iOM0yr-i,7;„„, = 12,000 K. 
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Fig. 7. — Comparison of observed and model Ha profiles for the 9 accretors with spectral 
types > M4. Note that the models do not include a wind component, and hence cannot 
reproduce the blueshifted absorption components seen in MHO-6 and IC348-415. Discrep- 
ancies near line center are likely due to rotation effects or a breakdown of the Sobolev 
approximation (see text). 
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Fig. 8. — Mass accretion rates as a function of stellar/substellar mass. Squares represent M 
values determined from Ha modeling and veiling in this work and Muzerolle et al. (2000a); 
solid circles are from Gullbring et al. (1998) (omitting binaries); pluses are from White & 
Ghez (2000) (including binary components); open circles are from White & Basri (2003). The 
arrows represent lower limits to M at 0.1 and 0.5 Mq measurable from Ha , as determined 
from the accretion models. 



